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We show that a 2 + 1 dimensional discrete surface growth model exhibiting KPZ class scaling can 
be mapped onto a two dimensional conserved lattice gas model of directed dimers. In case of KPZ 
height anisotropy the dimers follow driven diffusive motion. We confirm by numerical simulations 
that the scaling exponents of the dimer model are in agreement with those of the 2 + 1 dimensional 
KPZ class. This opens up the possibility of analyzing growth models via reaction-diffusion models, 
which allow much more efficient computer simulations. 
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I. INTRODUCTION 

The Kardar-Parisi-Zhang (KPZ) equation p] moti- 
vated by experimentally observed kinetic roughening has 
been the subject of large number of theoretical studies 
0, Q . Later it was found to model other important phys- 
ical phenomena such as randomly stirred fluid, Q , dissi- 
pative transport @,@], directed polymers and the mag- 
netic flux lines in superconductors It is a non-linear 
stochastic differential equation, which describes the dy- 
namics of growth processes in the thermodynamic limit 
specified by the height function /i(x, t) 

dth{x, t)=v + (tV^/i(x, t) + A(V/i(x, t)f + t]{x, t) . (1) 

Here v and A are the amplitudes of the mean and local 
growth velocity, cr is a smoothing surface tension coeffi- 
cient and 77 roughens the surface by a zero-average Gaus- 
sian noise field exhibiting the variance 



(77(x,t)77(x',i')> - 2D<5'^(x-x')(i-i') 



(2) 



Here d is used for the dimension of the surface, D for the 
noise amplitude and () denotes average over the noise dis- 
tribution. In 1 + 1 dimensions it is exactly solvable 0, 
but in higher dimensions only approximations are avail- 
able (see [9,] ) . In d > 1 spatial dimensions due to the com- 
petition of roughening and smoothing terms, models de- 
scribed by the KPZ equation exhibit a roughening phase 
transition between a weak-coupling regime (A < XO, gov- 
erned by A = Edwards- Wilkinson fixed point [lO|, and 
a strong coupling phase. The strong coupling fixed point 
is inaccessible by perturbative renormalization method. 
Therefore the KPZ phase space has been the subject of 
controversies and the value of the upper critical dimen- 
sion has been debated for a long time. Very recently non- 
perturbative renormalization and mode coupling theory 
has revealed a rich phase diagram, with more than one 
lines of fixed point solutions in the d - X space [111 ■ This 
suggests an upper critical dimension dc = A for KPZ, 



but an earlier numerical work |12| does not support this 
claim. 

In one dimension a discrete, restricted solid on solid 
realization of the KPZ growth is equivalent to the asym- 
metric simple exclusion process (ASEP) of particles 
[ll[l| (see Fig. [ID. 




FIG. 1: (Color online) Mapping of the 1-1-1 dimensional sur- 
face growth onto the Id ASEP model. Surface attachment 
(with probability p) and detachment (with probability q) cor- 
responds to anisotropic diffusion of particles (bullets) along 
the Id base space. 



In this discrete so-called 'roof-top' model the heights 
are quantized and the local derivatives can take the 
values Ah — ±1. By considering the up derivatives 
{Ah = 1) as particles and down ones as holes the rough- 
ening dynamics can be mapped onto a driven diffusive 
system of particles with single site occupancy. The ASEP 
model on the other hand is well known and its scaling 
properties are explored (for a recent review see p^). 

The extension of this kind of lattice-gas analogy to 
higher dimensions has not been considered to our knowl- 
edge. Instead hypercube stacking models were con- 
structed [H, [H, and surface configurations were 
mapped onto the d-state Potts spins defined on the sub- 
strate lattice itself. Especially 2-1-1 dimensional sur- 
faces were shown to be related to the six-vertex model 
with equal vertex energies [ll] and to the ground-state 
configurations of the anisotropic Ising model defined on 
the triangular lattice 19]. As a consequence the height- 
height correlation functions can be related to four-spin- 
correlation functions of the spin system. Very recently 



2 



the conformal invariance of the isoheight hnes has also 
been pointed out ^2Q] . 

Here we show that a 2 + 1 dimensional growth model 
exhibiting KPZ scaling can also be mapped onto a driven 
lattice gas. This is important from theoretical point of 
view, because the scaling behavior of driven diffusive sys- 
tem (DDS) has been studied intensively for a long time 
(for a review see ref. [2l|), thus results for DDS ma y b e 
exploited to understand KPZ better and vica versa [22] ■ 
Furthermore this conserved lattice gas and its generaliza- 
tions with anisotropics, disorder, or higher order terms 
can be studied effectively by bit-coded algorithms for ex- 
ample. 



II. MAPPING ONTO LATTICE GAS IN TWO 
DIMENSIONS 



As a generalization of the 1 -f 1 dimensional roof-top 
model, where the building blocks are squares let's put oc- 
tahedra on the square lattice, such that we get back the 
1-1-1 dimensional model in the x or y direction as shown 
on Figure [21 Surface adsorption or desorption events cor- 
respond to attachment or detachment of octahedra, re- 
spectively. The surface built up from the octahedra can 



the generalized Kawasaki updating rules 




FIG. 2: (Color online) Mapping of the 2-1-1 dimensional sur- 
face growth onto the 2d particle model (bullets). Surface at- 
tachment (with probability p) and detachment (with proba- 
bility q) corresponds to Kawasaki exchanges of particles, or 
to anisotropic diffusion of dimers in the bisectrix direction of 
the X and y axes. The crossing points of dashed lines show 
the base sub-lattice to be updated. Thick solid/dashed lines 
on the surface show the x/y cross-sections, corresponding to 
the Id model (Fig. [T]) 



be described by the edges meeting in the up/down mid- 
dle vertexes. The up edges in the x or y directions are 
represented by '-|-l'-s, while the down ones by '—1' in the 
model. In this way a single site deposition flips the four 
edges and means two '-1-1'^'— 1' (Kawasaki) exchanges: 
one in the x and one in the y direction. This can also be 
understood as a special 2d cellular automaton |23i| with 



-1 
-1 



1 -1 
1 -1 



(3) 



with probability p for attachment and probability q for 
detachment. We can also call the '-|-l'-s as particles and 
the '— 1-s as holes of the base square lattice. In this 
way an attachment/detachment update can be mapped 
onto a single step motion of an oriented dimer in the 
bisectrix direction of the x and y axes. To make a one-to- 
one mapping we update the neighborhood of sub-lattice 
points, which are denoted by the crossing-points of the 
dashed lines only. 

Since the three dimensional space can't be filled fully 
by octahedra, holes can occur among them, below the 
surface. Therefore this approximation of a surface growth 
may not sound to be faithful and the validity of KPZ 
growth rules requires confirmation. Note, however that 
in reality atoms are not cubes either and do not tile the 
three dimensional space completely. Furthermore very 
recently in bi-disperse ballistic deposition models [24, 25'!|, 
in which under-surface vacancies may occur KPZ scaling 
has been reported as well. 

The deterministic part of the KPZ equation ([1]), which 
can be obtained by averaging over the noise can be de- 
rived from the surface/dimer model similarly as it was 
done in 1 + 1 dimension [T3l|. If we apply the transfor- 
mation 

v(x,i) - V/i(x,t) (4) 

we get the Burgers equation for the height profile 

5tv(x, t) = crV2v(x, t) + Av(x, t)Vv(x, t) . (5) 

Our system is represented by two matrices A^, and Ay 
of sizes L X L, which contain discrete derivatives -1-1 or 
— 1 in a; and y direction, respectively (see Eqs. pU)) . ([TT|) ) . 
In two dimensions we introduce the vector variable Wij — 
{Ax{i~l,j),Ay{i,j—l)). This has the value (1, 1) in case 
of a dimer and (—1, —1) for a pair of holes. By setting 
up the master equation 

dtP{{a},t) = 5]^«:,,({a})P({a'},0 

for the probability distribution P{{'a},t), where the 
prime index denotes a state as a result of a generalized 
Kawasaki fiip Q the transition probability is given by 

— o P ^ o'j+ij + iCT.jj + A(CTi+1j + i — CTij) 



(1-^).- - ^2l 
(cr,j+ij+i X (T,,j) J 



(7) 



with A = 2^^^ — 1 parametrization. This formally looks 
like the one-dimensional Kawasaki exchange probability 
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(shown in jl3j), except the cross-product term, which is 
necessary to avoid surface discontinuity creation. The 
cross-product as a determinant cancels updates between 
configurations like (1, 1) —^ (1,-1). The nonlinear term 
vanishes for p = g (A = 0). The sign of the coefficient A 
of the nonlinear term can be interpreted as follows. For 
p > q positive nonlinearity (positive excess velocity) it is 
a consequence of growth with voids. 

To obtain Eq. ([5]) first one averages over the slope vec- 
tors 



(8) 



Then calculating its time derivative using the master 
equation the cross-product term drops out and one ob- 
tains 

+ A(CTij(CTi+lj + l — . (9) 

Here one can see the discrete second and first differentials 
ofWij corresponding to the operators of These differ- 
entials are one-dimensional because the dimer dynamics 
is also one-dimensional. Making a continuum limit in 
both directions and taking into account the relation of 
height and slope variables ^ we can arrive to the deter- 
ministic part of the KPZ equation ([T|). 

This agreement does not prove the equivalence of KPZ 
and the dimer model since they are just the first equa- 
tions in the hierarchy of equations for correlation func- 
tions. On the other hand from universal scaling point 
of view they show the equivalence of the leading order 
terms. We will show by numerical simulation that our 
mapping is faithful and reproduces the KPZ class sur- 
face growth behavior. 



III. THE SIMULATION ALGORITHM 

In the algorithm we extend the sequence of discrete 
slopes of the Id ASEP model (Fig. [1]) to local deriva- 
tives at sites in x and y directions of the surface 
(see Fig. [2]). The initially flat surface is presented as a 
regular sequence of '-|-l'-s and '— I'-s within both matri- 
ces. Periodic boundary conditions are applied to x and y 
direction. The system's evolution is simulated as follows. 

A site on the substrate plane is selected randomly. 
Then, we choose an attachment or detachment attempt 
according to their probabilities p and q. Generalized 
Kawasaki exchanges ^ of attachment or detachment are 
realized if 



or 



A,(i-l,j) A,(i,j) 
A,(i-l,j) A,{i,j) 

Ay{l,j-1) Ay(l,]) 



-1 1 
-1 1 

1 -1 
1 -1 



(10) 
, (11) 



respectively. Throughout this paper the time is mea- 
sured by Monte Carlo steps (MCS), i.e. L x L jump 



attempts correspond to one MCS. After certain time in- 
tervals data evaluation requires the reconstructions of the 
surface heights hx,y (t) by summing up the sequence of 
local slopes A^,, Ay. 



IV. RESULTS 

Starting from periodic, vertically striped particle dis- 
tribution, which corresponds to a flat initial surface we 
update the particle model by the rules deflned in the 
previous section. At certain time steps we calculate the 
hx.y{t) heights from the height differences Ax.y The 
morphology of a growing surface is usually characterized 
by its width 



WiL,t) 
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1/2 



(12) 

In the absence of any characteristic length, growth pro- 
cesses are expected to show power-law behavior of the 
correlation functions in space and height and the surface 
is described by the Family- Vicsek scaling ^26(| 



W{L,t) oc tf^, for to«t«ts 
oc L°, for t » ts . 



(13) 
(14) 



Here a is the roughness exponent and characterizes the 
deviation from a flat surface in the stationary regime 
{t >> tg), in which the correlation length has exceeded 
the linear system size L; and (3 is the surface growth 
exponent, which describes the time evolution for earlier 
(non-microscopic t » to) times. The dynamical expo- 
nent z can be expressed by the ratio 



z = a/P 



(15) 



In case of up-down symmetry (p = 1, g = 1) the non- 
linear term is dropped, and the KPZ equation ffl simpli- 
fies to the Edwards- Wilkinson (EW) equation [13] ■ Since 
the upper critical dimension of this equation is: dc = 2, 
mean-field behavior, characterized by a = /3 = and log- 
arithmic scaling is expected by field theory. Indeed, the 
width of the surface grows like 



W^it) = aln(t) +6 



(16) 



as shown in Fig. [31 The prefactor a obtained by fitting 
the L 1024 curve in the 20 < t < 1000 region with the 
form is a = 0.152(8). This is in agreement with the 
theoretical estimate for the EW equation D/iAira) [27] if 
take into account the exact value for the stiffness constant 
(or surface tension): a/D = 7r/9. This constant was 
identified by [l^ through the correspondence between 
the exact calculation of the four-spin correlation function 
of the zero-temperature triangular Ising antiferromagnet 
[2^ and the discrete height-height correlation function in 
real space in the interface model. A factor 5- = 2/3cr is 
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FIG. 3: (Color online) Logarithmic surface growth in 
case of up-down symmetry for different sizes L — 
64, 128, 256, 512, 1024 (bottom to top). The dashed line shows 
the fitting with the form (|16p . Inset: width saturation values 
for different system sizes L in the long time limit. 



of '30'|. Using these surface exponents and the scaling 
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FIG. 4: (Color online) Scaling collapse for p — 1, q — with 
2 + 1 dimensional KPZ class exponents for different sizes: 
L = 64, 128, 256, 512, 1024 (bottom to top). 



coming from the 2/-\/3 triangular lattice site per surface 
element and the l/\/3 of the octahedron/cube surface 
fraction, thus the theoretical estimate is: a ~ 0.151981. 

The saturation values are expected to exhibit logarith- 
mic growth 

W^iinf.) = lim W^{t) = cln(L) + d (17) 

t — ^oo 

with the system size [13]. As can be seen in the inset of 
Fig. [3]this really happens with the prefactor c = 0.30(1), 
which agrees with the theoretical value c — 2a ~ 0.304 
again. 

For pure deposition p = 1, g = 0, or in case of other 
general up/down asymmetric cases, we saw power-law in- 
crease of the surface width, in agreement with the scaling 
hypothesis (see Fig. 2]). For the the largest system 
that we have investigated (L = 1024) we fitted W{t) in 
the 100 < t < 10000 time window with a power-law and 
obtained f3 — 0.23(1). This value agrees quite well with 
the numerical estimate for the 2 -|- 1 dimensional KPZ 
class (/3 = 0.24) provided in ref. 0]. Note however that 
for smaller system sizes the exponent estimate is some- 
what smaller, due to corrections to scaling, but one can 
clearly see a convergence towards higher values and a 
better collapse as L ^ oo. Large scale simulations with 
an effective, bit-coded version of our algorithm could re- 
sult in very precise estimates. The systematic tendency 
towards an asymptotic behavior has been found in Fig.|4l 

The saturation values W{inf.) for different system 
sizes also scale well with (ITil) and with the exponent 
a = 0.38(1) of the 2+1 dimensional KPZ class As- 
suming corrections to scaling of the form W ~ A2L°'{1 + 
B2L~'^ the fitting to our data resulted in very small ef- 
fect: a = 0.377(15), which marginally overlaps with the 
value of [1^ but does not support the proposal a = 2/5 



law (jlSp we estimated the dynamical exponent to be: 
z = 1.64(1), which is somewhat greater than what one 
finds for the 2 + 1 dimensional KPZ class in Q (z = 1.58). 
We think that this is due to the correction to scaling ob- 
served in the time dependence discussed above. If we 
scale the time with the dynamical exponent z = 1.64 
we obtain a good scaling collapse of the growth data for 
different sizes (FigH]) in agreement with the p3ll4p law 
again. Our exponent estimates also satisfy the a + z ~ 2 
scaling relation within error margin. This implies that 
the Galilean invariance holds and the lattice model in- 
deed lies in the 2 + 1 dimensional KPZ universality class. 

V. CONCLUSIONS AND OUTLOOK 

We have pointed out the possibility of mapping of a dis- 
crete surface growth processes onto a conserved, driven 
lattice gas model of oriented dimcrs, which move perpen- 
dicularly in two dimensions. The straight line motion of 
dimers in the two dimensional space is very similar to the 
motion of particles of the ASEP process. The difference 
is that since the dimers are extended objects, their mo- 
tion is slowed down by the dimer particle exclusion and 
the sub-lattice update as compared to the single particles 
of the ASEP. As a consequence their motion is described 
by somewhat larger dynamical exponent {z ~ 1.64) than 
that of the ASEP {z = 3/2), so the change of z{d) seems 
to be a purely topological phenomena in KPZ. This pro- 
vides a better understanding of the relation of univer- 
sality classes of surface classes to those of the reaction- 
diffusion models, m, [111 . Interestingly the cc/y sym- 
metric surface dynamics maps onto a strongly anisotropic 
reaction-diffusion model. 

We have found KPZ or EW scaling by numerical sim- 
ulations, hence we showed that lattice anisotropy and 
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under-surface vacancies are irrelevant. Our simulation 
results for the 2 + 1 dimensional EW case reproduced the 
theoretically expected logarithmic scaling, with the cor- 
rect leading order coefficients. For the KPZ scaling our 
roughness exponent result is in the middle of the range 
obtained by various numerical exponent estimates: i.e. 
between a = 0.36 [s^, and the field theoretical value 
a = 0.4 d^]. Our a = 0.377(15) coincides with that of 
the numerical study 13611 and agrees with the renormaliza- 
tion results a = 0.38 [33|. It overlaps marginally with the 
simulation results a = 0.393(3) [20| as well. Our growth 
exponent estimate /3 = 0.23(1) matches the results of [13] 
(/3 = 0.221(2)) and (/3 = 0.229(5)), obtained by in- 
dependent numerical fitting procedures. The dynamical 
exponent of this study is also in the range provided in 

Our model provides an efficient way of simulations and 
opens us the possibility to study more complex growth 
models relevant in recent interest of self-organizing sur- 
face nanosystem An optimized, bit-coded version 



of our code, which manipulates the two-dimensional bit- 
field by logical operations runs roughly 10 times faster 
than the current version and will be published elsewhere. 
For example the Bradley-Harper 39] and the debated 
Kuramoto-Sivashinsky [40] models with their modifica- 
tions can be investigated numerically and are the subject 
of a forthcoming publications. 
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